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Surface modes are a special type of wave localized at the interface separating 
two different media. In periodic systems, the modes localized at the surfaces 
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are known as Tamm states [1], first found as localized electronic states at the 
edge of a truncated periodic potential. Surface states have been studied in 
different fields of physics, including optics [2], where such waves are confined 
to the interface between periodic and homogeneous media. 

Photonic crystals, artificially fabricated periodic structures with bandgap spec- 
tra [3], can be used for controlling the properties of light in different devices 
including dielectric mirrors and waveguides. In many such applications, pho- 
tonic crystals are finite, and are terminated at surfaces where electromagnetic 
waves are significantly affected by the breaking of the translational invariance 
by the underlying periodic structure. 

Although the properties of electromagnetic waves in bulk photonic crystals are 
well understood, studies of electromagnetic waves near the surfaces of photonic 
crystals are relatively limited and, unlike electronic Tamm states, a truncated 
photonic crystal does not always support surface states [4]. Surface states in 
photonic crystals have only been shown to exist under appropriate changes to 
the surface layer, such as a termination though the surface cell or a change to 
the surface geometry or material properties [5,6,7]. Furthermore, such surface 
waves are known to be significantly sensitive to the surface termination [6,8]. 

In this Communication, we study surface waves in photonic crystals with ter- 
minated (but not altered) surface structures. For the first time to our knowl- 
edge, we show the existence of strongly localized surface modes in binary 
photonic crystals, in a complete analogy to Tamm states in electronic sys- 
tems. Our analysis of the linear properties of the diatomic structure illus- 
trates the dispersion relations and field localization of the surface states, and 
the subsequent influence of variations to the crystal geometry. Additionally, 



we study the nonlinear behavior of the surface states when the dielectric func- 
tion of the photonic crystal includes a Kerr style nonlinearity, and discuss the 
nonlinearity-induced tunability of nonlinear Tamm surface states. 



2 Model and numerical methods 

We consider the propagation of the TM-polarized waves in a two-dimensional 
binary photonic crystal formed by two square lattices of larger and smaller 
dielectric rods, as shown in Fig. 1(a), where r^ and r^ are the rod radii popu- 
lating the lattices. These two elemental lattices merge to form a single lattice 
of period a, with a unit cell created by a single, large rod at the center of the 
cell and a single small rod distributed as four quadrants at the corners of the 
cell; thereby forming a unit cell with two rods. In general, the binary structure 
reduces the geometric symmetries of the photonic crystal, thereby lifting the 
degeneracies of some of the crystal states [9]. However, the structure we study 
here does not possess this capacity and as such maintains the C^^ point group 
symmetry of the constituent square lattices. 

We consider the surface termination that has translational symmetry along 
the surface where a row of either large radius rods r^ or small radius rods r?, 
form the surface layer exclusively. In all cases discussed below only complete 
rods are considered within the surface layers, similar to the original context 
of Tamm states. 

For these binary crystals, we set the radius of the rods as r^ = 0.21a and 
Tf, = 0.13a. The dielectric rods of the photonic crystal correspond to high 
potential optical regions. In this study, we set the linear dielectric constant of 
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Fig. 1. (a) Structure of a binary photonic crystal with a surface termination. Shad- 
ings (red and green) show the equivalent irreducible Brillouin zone of an infinite 
and semi-infinite photonic crystal projected onto the crystal geometry, (b) Spatial 
profile of the linear surface mode as the color plot of l-B^P- 
the rods to Sr = 11.56, corresponding to AlGaAs at a wavelength of 1550nm, 

and assume negligible losses. 

To analyze the surface states, we employ two complementary numerical meth- 
ods: the plane wave expansion method (PWE) [10,11] and the finite-difference 
time-domain method (FDTD) [12,13]. Due to the normal computational in- 
tensive nature of the FDTD method and the need to reduce the time step 
to maintain numerical stability [14] for the nonlinear case, the method is ex- 
tremely time consuming. To reduce this burden, we model a semi-infinite crys- 
tal using a supercell representation [15,13] that has a transverse size of one 
unit cell, and contains nine unit cells of the photonic crystal, including the sur- 
face cell, and nine unit cells of vacuum in the lateral direction. The supercell 
is bordered by a perfectly matched layer [16]. The boundaries perpendicular 
to the surface are periodic, and configured with a vr spatial phase shift that 
sets the wave vector to the edge of the surface Brillouin zone. 

The supercell is excited using a two-dimensional Gaussian pulse, centrally 
located within the spatial domain and the photonic band gap, that energizes 
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Fig. 2. (a) Surface mode dispersion relationships for the surface formed by larger 

(upper solid line) and smaller (lower dashed line) rods. (b,c) Photonic band gap 

and surface mode frequencies with respect to dielectric rod radii. In (b) the larger 

rods of radius r^ are held constant while the radius r^, of the smaller rods is varied, 

whereas in (c) the smaller rods of radius r^ are held constant while the radius r^ of 

the larger rods is varied. 

the resonant modes of the crystal and the surface structures. A discrete Fourier 

transform of the time evolution data, recorded at the surface of the photonic 

crystal, is used to determine the spectral response of the surface states. In a 

similar manner, the spatial mode profile is calculated from the time evolution 

data of the full simulation domain at the peak frequency of the surface state. 

3 Linear Tamm states 



For the linear surface states, the dispersion relationships are presented in 
Fig. 2(a). The dark grey region within the dispersion diagrams depicts the 
continuum of free-space states that exist above the line light {uo = ck) that 
are not bound to the surface. Light grey regions within the diagram indicate 
the Bloch states of the infinite photonic crystal; states within these regions 



can couple into the crystal, and again are not bound to the surface. 

Surface states are formed within the second photonic bandgap of the infinite 
crystal, and are shown with a solid black line in the diagram. The surface 
states occurs within the second band gap when the surface is terminated in 
a row of large rods of radius r^. A localized surface state can form within 
the first bandgap when the crystal is terminated in a row of smaller rods, as 
indicated by the dashed line in Fig. 2(a). This state, however, can couple to 
the Bloch states of the infinite crystal, and we do not consider it in detail. 

We analyze the characteristics of the surface states by varying the radii of the 
constituent rods of the photonic crystal. Figures 2(b,c) show the influence of 
the rod radii on the bandgap edges and surface state frequencies. The larger 
rods principally define the upper bandgap edge for the T — X direction and 
have virtually no effect on the lower bandgap edge, whereas the smaller rods 
set the lower bandgap edge with very weak effect on the upper bandgap edge. 
Figures 2(b,c) also confirm that the surface mode does not exist when the 
radii of the two constituent rods are equal 

The consistency of the location of the surface states within the bandgaps 
signifies the states existence domain is predominantly defined by the crystal 
geometry; a feature attributed to the binary nature of the crystal, which has 
previously been shown to provide robust spectral features in the presence of 
geometric variations [9]. The spatial field localizations can be understood by 
noting the symmetries of the Bloch waves above and below the bandgap. In 
the r — X direction below the second bandgap, the Bloch wave symmetries 
place the electric field within the smaller rods, whereas for above the bandgap 
the field resides between the rods in free-space, with some overlap into the 
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larger rods. As the surface is formed by large rods and the states are near 
the upper bandgap edge, they take on the Bloch mode symmetries of this 
frequency region. 

4 Nonlinear Tamm states 

We study the properties of nonlinear surface states through the addition 
of a third-order susceptibility term within the polarization field of all rods. 
Throughout our analysis we only modify the x''^^ coefficient and maintain the 
input field intensity and density. Under steady-state conditions, the principal 
effect of the nonlinear term is to induce an intensity dependent change to the 
dielectric strength of the surface rods. Simulations of the nonlinear dielectric 
rods of the photonic crystals are performed using the FDTD method [17,12]. 
This is achieved by adding a x^^^ term to the polarization as: P = x'-^'^-'E -|- 
;\;(^)|EpE. The analysis of the linear surface states revealed that the field in- 
tensity varies appreciably across the surface rods. For the nonlinear surface 
rods this causes an equivalent change in the dielectric constant. In our anal- 
ysis of the nonlinear surface states we consider only a focussing nonlinearity 
( X^^^ > 0) that results in localized increases in the dielectric index which are 
proportional to the cube of the electric field. In turn, this leads to a decrease 
in the resonant frequency of the surface state. 

Figure 3(a) illustrates the surface mode frequency shift as the strength of 
the third-order susceptibility grows. A strong frequency shift is observed for 
moderate changes in x^\ with a threshold to the onset of strong frequency 
shifting at approximately x''^^ = 1 x lO^^fim^/V. These x*"^"* values exemplify 
the significant dynamic range of the nonlinear surface state of a binary pho- 
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Fig. 3. (a) Nonlinear frequency shift of the surface state vs. the nonhnear suscep- 
tibihty for r^ = 0.21a and r^ = 0.13a; (b,c) Spatial profile of the surface modes 
as the color plots of the field density |i?2p for x = 2 x 10^"^ p.m'^ /V (left), and 
;^;(3) = 9 X IQ-'^pm^/V (right). 

tonic crystal using realistic nonlinearities; for example a typical value of x^^^ 
for AlGaAs is x^^^ = 8.2 x IQ-^ fim^/V [18]. 

The distinct spatial field profile of the surface state, and its dramatic change 
with dielectric strength leads to the strong x^^^ sensitivity. This is demon- 
strated in Fig. 3(b,c), which shows the spatial field profiles for (b) x^^^ = 
2 X lO^^nm^ /V, and (c) x^'^^ = 9 x IQ-^ixm^ /V. In a homogeneous nonlinear 
focusing medium, a nonlinearity-induced change in the dielectric constant re- 
sults in a greater localization of the light, which under appropriate conditions 
can lead to spatial surface solitons. However, for the surface state of the bi- 
nary photonic crystal, the increased dielectric constant causes the conditions 
for the surface resonance to diminish, resulting in a reduction in the field 



localization, leading to a complex balance between the nonlinearity and sur- 
face interaction. In general, as the nonlinear surface state evolves and reaches 
steady state conditions, the field intensity within the surface rods reduces. This 
effect provides optical limiting that prevents saturation of the nonlinearity and 
unrealistic changes to the dielectric strength. Additionally, the competing ef- 
fects due to nonlinearity and the surface can lead to complex behavior, as seen 
for x''^^ = 7 X 10^^ fim'^/V where the frequency of the surface state increases. 
Another consequence of the nonlinear surface structure is the formation of a 
near-surface defect state, as seen in Fig. 3(c), where the localized state forms 
in the neighboring row of small rods. 

Recent studies have demonstrated the application of optical surface states of 
photonic crystals for guiding light, the formation of high-quality micro cavities, 
and sub-wavelength imaging [19,20,21]. We expect that the nonlinear tunable 
Tamm states found here would provide substantial nonlinearity-induced con- 
trol within these applications leading to novel all-optical surface devices such 
as optical limiters and switches. 



5 Conclusion 

In conclusion, we have predicted the existence of surface Tamm states in bi- 
nary two-dimensional photonic crystals. Using nonlinear surface rods, we have 
highlighted the dynamic tunability of the surface states resulting from the 
unusual spatial distribution of the surface mode, and its geometric transfor- 
mation through a competition of the nonlinearity and surface effects. 
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